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. Abstract 



We consider matrix-model representations of the meander problem which de- 
scribes, in particular, combinatorics for foldings of closed polymer chains. We intro- 
duce a supersymmetric matrix model for describing the principal meander numbers. 
This model is of the type proposed by Marinari and Parisi for discretizing a super- 
string in D = 1 while the supersymmetry is realized in L> = as a rotational sym- 
metry between bosonic and fermionic matrices. Using non-commutative sources, we 
reformulate the meander problem in a Boltzmannian Fock space whose annihilation 
and creation operators obey the Cuntz algebra. We discuss also the relation be- 
tween the matrix models describing the meander problem and the Kazakov-Migdal 
model on a I?-dimensional lattice. 
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1 Introduction 



Matrix models is a standard tool for describing discretized random surfaces (or, equiva- 
lently, strings) [Q]. A supersymmetric extension of this construction was first proposed by 
Marinari and Parisi [0] and studied for the D = 1 dimensional target space. 

We introduce in the present paper supersymmetric matrix models in the D = di- 
mensional target space which differ from the Hermitean supermatrix models of Ref. 0] 
and are a version of the Marinari-Parisi construction in D = 0. Our model deals with 
the "superfields" 

Wa={B,F), Wa={B\F) (1.1) 

where a = 1, 2 while B and F are complex bosonic and fermionic {i.e. Grassmann valued) 
N X N matrices, respectively. 

Since the propagators for both bosonic and fermionic matrices coincide: 

the supersymmetry reduces in D = simply to rotations between the B- and F-compo- 
nents. The proper transformation reads 

6B = eF, 6F= -eB , 

6B^ = Fe, 6F = -Bh, (1.3) 

where e and e are Grassmann valued. 

Any potential, which is symmetrically constructed from the "superfields" (|1.1|) , is 
supersymmetric so that contributions from the loops of the bosonic and fermionic matrix 
fields are mutually cancelled which is the key property of the supersymmetry. The simplest 
Gaussian supersymmetric potential reads 

2 _ 

Kiauss = NJ2 trWaWa = N ti [b^ B + F f) , (1.4) 

a=l 

which reproduces the propagators ( |1.2| ). It is obviously invariant under the rotation (|1.3|) 
even when e and e are fermionic N x N matrices. It is also clear from Eq. ( |1.4| ) why one 
needs complex matrices in D = 0: the trace of the square of a fermionic matrix vanishes. 

We elaborate in this paper the technique for dealing with the supersymmetric matrix 
models on an example of the one which describes combinatorics of the meander num- 
bers. This challenging combinatorial problem, which is described in Sect. 2, is not yet 
solved. The associated matrix model, which describes a physical problem of enumerat- 
ing different ways of foldings of a closed polymer chain, is of a next level of complexity 
with respect to the Hermitean one- or two-matrix models, the multi-matrix chain and its 
multi-dimensional extension — the Kazakov-Migdal model [Q. 
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We introduce in Sect. 2 the complex matrix model which is equivalent to the Hermitean 
one in describing the meander numbers. We construct then the supersymmetric matrix 
model which describes the principal meanders. We discuss also the relation between the 
matrix models describing the meander problem and the Kazakov-Migdal model on a 
D-dimensional lattice. 

Using non-commutative sources, we reformulate in Sect. 3 the meander problem as 
a problem of averaging in a Boltzmannian Fock space whose annihilation and creation 
operators obey the Cuntz algebra. The averaging expression is represented in the form of 
a product of two continued fractions. 

The Appendix A contains a solution of the combinatorial problem of summing over 
words built up of unitary matrices, which is equivalent to the Kazakov-Migdal model 
with the Gaussian potential, via free random variables. 

In the Appendix B we demonstrate how the equations of Sect. 3, which are obtained 
using the matrix-model representations, can be alternatively derived pure combinatorially. 

We comment in the Appendix C on a possibility of solving the meander problem 
via free random variables. We show that this approach does not work for the meander 
problem since the variables are not free for this case so that the theorem of addition of 
free random variables is not applicable. 



2 Matrix models for the meander problem 

The meander problem is known to people working on Quantum Field Theory since the 
middle of the eighties from V. Arnold. The problem is to calculate combinatorial numbers 
associated with the crossings of an infinite river (Meander) and a closed road by 2n 
bridges. Neither the river nor the road intersects with itself. These principle meander 
numbers, M„, obviously describe the number of different foldings of a closed strip of 2n 
stamps or of a closed polymer chain. 

One can consider also a generalized problem of the multi-component meander numbers 
M^'^) which are associated with k closed loops of the road so that M„ = M^^^ . The results 
of a computer enumeration of the meander numbers are presented in Refs. ^ up to 
n = 12. 

2.1 Hermitean matrix model for meanders 

Meanders can be described by the following Hermitean matrix model 

^See Ref. ||^ for an introduction to the subject. 
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where the integration goes over the N x N Hermitean matrices Wa {a = 1, . . . ,m) and 

0. The logarithm in Eq. leaves only one closed loop of the field 0. The coupling 
constant c is associated with the (quartic) interaction between Wa and 0. 

Expanding the generating function ( |2.1| ) in c and identifying the diagrams with the 
ones for the meanders, one relates the large- limit of J-'nxn{c) with the following sum 
over the meander numbers 

oo 2n n 

hm :F^Mc) = T.i^Y. Mi'^^' . (2.2) 

n=l 2^ k=l 

The N ^ oo limit is needed to keep only planar diagrams as in the meander problem. 

The RHS of Eq. (^]l|) can be expressed entirely via the Gaussian averages of W^s. 
This leads to the following representation of the meander numbers: 

^MWm'== {^^^Wa,Wa,---Wa,^.,Wa,S)^ , (2-3) 

k=l ai,a2,---,a2n-i,a2n=l ^ / GauSS 

where the average over W^s is calculated with the Gaussian weight — the same as in ( |2.1| ). 
This formula can be proven by calculating the Gaussian integral over in Eq. ( p.l| ), 
expanding the result in c and comparing with the RHS of Eq. ( p. 21 ) . The factorization at 
large is also used. 

The principle meander numbers M„ are given by Eq. ( |2.3| ) as the linear-in-m-terms, 

1. e. as linear terms of the expansion in m. This looks like the replica trick which suppresses 
higher loops of the field W . 

The ordered but cyclic-symmetric sequence of indices ai, 02, . . . , a2n-i, a2n is often 
called a word constructed of m letters. The average on the RHS of Eq. (p.3|) is the 
meaning of a word. Thus, the meander problem is equivalent to summing the squares of 
all the words with the Gaussian meaning. 

The Gaussian averages on the RHS of Eq. ( |2.2| ) can be represented, making the Wick 
pairing, via the Kronecker deltas: 

/ 1 



a.2n-ia2n ~^ plauar permutations . 
\-'V / Gauss 

(2.4) 

The "planar permutations" means here that one should sum up over all the permutations 
of the indices Oj's which are consistent with the planarity. This is standard for the large- iV 
limit. 

To calculate the meander numbers, one should sum up the square of the RHS of 
Eq. ( p.4| ) over Oj's as is prescribed by Eq. ( |2.3| ). This is a convenient practical way of 
calculating the meander numbers. 

Since for m = 1 

\n /Gauss {n + l)\n\ ^ ' 
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which is known as the Catalan number of the order n, one gets from Eq. (p.3|) 

n 

E Ml^^ = Cl . (2.6) 

k=l 

This is nothing but the first sum rule of Ref. 

2.2 Complex matrix model for meanders 



The combination of deltas on the RHS of Eq. (|2.4| ) can be alternatively represented as 
the Gaussian average over the complex matrices: 

= y n dWldWae-^^:-^ '-'^^'^^j.i^W^M^ ■ ■ ■l^a.„_,W^i„ . (2.7) 

a=l 

The generating function, associated with the representation of the meanders via the 
complex matrices, reads 



N'^\\J ) I Gauss 

Ar2 



^ /ndIVtrfiy.e-^i::^i-^^^^Mn(| d0id</>2e-^) (2.8) 



with 

Af JV ™ 

^ = -tr02+ tr ciV ^tr (^iW^.Wa) • (2.9) 

a=l 

Here 01 and 02 are Hermitean while Wa (a = 1, . . . , m) are general complex matrices. 

Quite similarly to Eq. ( p.l|) where the Hermitean matrix can be represented in a 
diagonal form = diag the matrices 0i and 02 in Eqs. (|2.8|) , (|2.9|) can 

always be made diagonal: 

01 = Ai^diag (aW,...,aS^)), 

02 = A2 = diag(A«,...,Af^). (2.10) 
This can be shown representing 0i and 02 as 

01 = l^lAi^^i , 02 = ^]^A2^^2 (2.11) 
and absorbing the unitary matrices Qi and Q2 by the transformation of Wa'- 

Wa nlWafii , — > n\w^n2 ■ (2.12) 

The measure dW^dWa does not change under the transformation ( p.l2| ) since Wa are 
general complex matrices. 
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It is convenient to introduce one more generating function 

where only one component of Wa, say the first one, enters the averaging expression. 
Differentiating the generating function ( |2.8D with respect to c and noting that all m 
components of Wa are on equal footing, we get the relation 

c^^=mM{c) (2.14) 

between the two generating functions. 

In order to show how the complex matrix model recovers the meander numbers, let us 
replace 0*/ or 02"' iii the numerator of Eq. (|2.13 ) by N~^d/d(l){^ or N~^d / d(f)2 ■, respectively. 



and integrate by parts. Repeating this procedure iteratively, we get 

oo n 

M{c) = Y.c'-Y.Mi'^^'-\ (2.15) 

n=l k=l 

with 

n m , -I 1 2 

j:Mi'W-'= E {^i^W,Wl---Wa,„,,Wl)^ , (2.16) 

k=l a2,-..,a2„_i,a2„=l ^ / GauSS 

where the Gaussian average is defined by Eq. (12.7]). Equation ( |2.16| ) can be alternatively 
derived calculating the Gaussian integrals over (pi and 02 in Eq. ( 2.13|) by virtue of 



M02e-^ = /rf02e-^*'^*^+^'='^S--^^'^(*^^''^^'*^^'^") 



det-i/2 



m ^ 

I ® I - E WaW^ ® (WkW^) 



a,b=l 



[2.17) 



(^.tT{WWY) =C„, (2.18) 

\jV / Gauss 



For m = 1 the formula 



which is analogous to Eq. ( p.5|) , holds for the complex matrices. This results again in 
Eq. (IH). 

It is instructive to consider also the case when W is a. fermionic Grassmann valued 
matrix d la Ref. We shall denote the fermionic matrix as F and its conjugate as F. 
Then we get 

ItrfFFrX -/ ^ = 2p(even) 

N^'y^^J /Gauss ~\C, n = 2p+l (odd) " ^'''^^ 



^See Ref. M for a review 



6 



Since each loop of the fermionic field is accompanied by a factor of (—1), we arrive at the 
sum rule 

This is nothing but the second sum rule of Ref . . 

Note that the trace of the square of a fermionic matrix vanishes because of the anticom- 
mutation relation imposed on the components. This is why we did not consider Hermitean 
fermionic matrices and used first a representation of meanders in terms of complex ma- 
trices to discuss fermionic representation of meanders. Fermionic matrix models are a 
natural representation of the notion of the signature of arch configurations of Ref. M. 



2.3 General matrix model for meanders 

The consideration of the previous Subsection suggests the following representation of the 
meander numbers via a general complex matrix model which includes both bosonic and 
fermionic matrices. 

Let us consider a general Wa which involves both bosonic (complex) and fermionic 
(Grassmann) components: 

Wa = ( Bi, B2, . . . , Bm^, Fi, F2, . . . , Fmj, j , 



Wa = (Bl Bl,..., Bl. Fi, A, . . . , J . (2.21 



Here and rrif {m = mh+mf) are the numbers of the bosonic and fermionic components, 
respectively. Let us define the generating function J^{c) by the formulas 



Gauss 



and 



= J^J U dWadWa e -^Sr^i ^■-^''^'^ In (^J d^cp^ e j (2.22) 



S = -tT<Pl + -tT<Pl-cNY.tT {<PlWa<p2Wa) , (2.23) 

a=l 



which generalize Eqs. ( p^.8| ) and ( |2.9| ). Then, the generating function ( |2.22| ) is related to 
the meander numbers by 



oo ^2n n 



^hn ^^x^(c) = E ^ E M^n\^^ - ^ff ■ (2-24) 



fc=i 



Here mj emerges with the minus sign since fermion loops are always accompanied with 
the minus sign. 



7 



Equation ( |2.16|) is extended to the given general case of both bosonic and fermionic 
matrices as 



a2,---,a2n-i,a.2n=i 



k=l 



1 1 



\N"^"""'' " "^"-^ ' / Gauss \iV 
Analogously, Eqs. (|2.13|) , (|2.14|) and ( |2.15|) are generalized as 



(lt^Wa,„Wa,„_,---Wa,W,) . (2.25) 

Gauss 



M(c) = c ' 



Gauss 



"^"^c"^ = {rrib - mf)M{c) 



oo n 



(2.26) 
(2.27) 

(2.28) 



where 

M(c) = ^ c^" ^ Mi'\mh - rriff-^ 

n=l k=l 

Equations ( |2.13|) , ( |2.14 ) and ( |2.15 ) are obviously reproduced when rrif = 0. 

In order to prove Eqs. ( |2.25| ), ( |2.28| ), let us first note that Eq. ( p. 17] ) is extended to 
the general case of both bosonic and fermionic components as 



hd(P2 e 



92 e 



92 e 



det-i/2 



where 



aia) 



-Y^'^'t''2 + h''^T,Zb=l {Wa<p2WaWb4>2Wb) 
-f tr <t>l + lc^NY,Zb=l '^(") *^ i^2WaWt<l>2WtWa) 



a,b=l 

f 1 for 5 
I -1 for F 



[2.29) 



[2.30) 



is the signature factor of the component Wa. 

Expanding the determinant ( p.29| ) on the RHS of Eq. ( p.22| ) in c^, we get the following 
representation 

oo ^2n m 

lim J^nxn{c) = — Y] cr{ai)a{a3) ■ ■■a{a2n-i) 

w— >oo zn 1 

n=l ai,---,a2„-i,a2„=l 

\jV /Gauss \jV /Gauss 

oo 2n m 

= V— V 

^ 2n ^ 

n=l ai,---,a2n-i>a2n=l 
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where the signs of the fermionic components are transformed at the last step using 

a(ai)a(a3) ■ ■ ■cr(a2n-i) tr Wa.Wa^ ■ ■ ■ Wa,„_,Wa,,}j ^^^^ 



.j^tTWa,Wa,---Wa,„_,Wa,) • (2.32) 



Therefore, the order of matrices in Eq. ( 2.25 ) is chosen in the way to absorb the signature 
factors for the fermionic components. 

2.4 Supersymmetric matrix model for principle meander 



Having the representation ( |2.25|) of meanders via general complex matrices (either bosonic 



or fermionic), we can utilize the idea of supersymmetry to kill the loops of the VT-field 
instead of the replica trick. Let us consider the two-component Wa whose first component 
is bosonic while the second one is a fermionic matrix: 



Wa = {B, F) , Wa^ Wl = [B\ F) , (2.33) 
that is given by Eq. ( |2.21|) with = rrif = 1. 



The generating function ( p.22| ) equals zero for the supersymmetric model since all 
the loops of the B and F fields are mutually cancelled. One should use alternatively the 
generating function ( |2.26| ) which can be represented for the supersymmetric matrix model 



M{c) = (^tiBBUnf f d(j)id(j)2e-^]) . (2.34) 
\N \J } I Gauss 

Here S is explicitly given by 



^ = ^ tr 02 + ^ tr 0^ _ cA^ tr (015^2^) - cN tr (0iF02F) (2.35) 



as is prescribed by Eq. ( |2.23| ) with Wa substituted according to Eq. ( |2.33| ). The equivalence 
of Eqs. ( |2.26| ) and ( [^.34| ) in the supersymmetric case can be proven replacing B in the 
integrand on the RHS of Eq. ( |2.34| ) by N~^d / dB\ integrating by parts, and recalling that 



J^(c) = 0. 

All the multi-component meanders in Eqs. ( p.25|) , (|2.28|) vanish in the supersymmetric 
case and we get the following representation for the principle meander 

oo 

M{c) = c'^Mn (2.36) 

n=l 

with 

2 

a2,---,a2n-l,Cl2n = i 

X (^trBWa,---Wa,^_,Wa,\ ( ^ tr 1^,,^, ■ ■ ■ 1^,, ,(2.37) 
\ -/v / Gauss \ -iV / Gauss 
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n 


Structure 


Combinatorics 


Value 


Contribution to M„ 


1 


\ 1 


1 


1 


1 


2 


m 

{B^F^) 


1 ] 
2 




2 
1 


4 1 


2 






(B^) 


1 ' 




5 


25 ' 








{B^F'^) 


4 




2 


— 16 






3 


(B^FB^F) 
(BF'^BF^) 


2 



1 


> ID 


1 

1 


_2 

n 
1 








\^ 1 


1 ' 




14 


196 








(B^F^) 


6 




5 


-150 








(B^FB^F) 


6 




2 


-24 








(B^F^) 


4 














4 


(B^F^B^F^) 


2 


>35 


1 


2 




► 42 


(B^F'^BF^) 

{B^FB^F^) 

{B^F^) 

(BF^BF^) 

(BF^BFB^F) 


4 
4 
2 
2 
4 


2 

1 

1 


16 


-2 



4 





Table 1: Calculation of the principle meander numbers in the supersymmetric matrix 
model ( |2.37] ) up to n = 4. 



where we kept trace of the order of matrices how it appears from Eq. (|2.25|) . The 
signs, which are essential for the fermionic components, have been transformed using 
the formula (|2.32 ) with a [a) being the signature factor of the component Wa defined by 



Eq. (|30D. 

Equation ( |2.37| ) is a nice representation of the principle meander numbers which looks 
more natural than the one based on the replica trick. A hope is that it will be simpler 
to solve the m = 2 supersymmetric model than a pure bosonic one at arbitrary m. How 
the representation ( |2.37] ) reproduces the principle meander numbers is illustrated by the 
Table |I| up to n = 4. 

Alternatively, one can calculate 

a2,---,l2n-l,a2n = l 

1 t.F«'.. . . . W„„_.W„,)^_ (1 tHy„„,y„„_. . . . . (2.38) 

The results are presented in the Table |^. 

The analogous results for the pure bosonic case given by Eq. ( p.l6|) at m = 2 are 
presented in the Table 3- 
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n 


Structure 


Combinatorics 


Value 


Contribution to — M„ 


1 


\ / 


1 


1 


-1 






1 ] 







1 




2 


\ 1 


2 




1 


-2 J 








1 ' 




1 


-1 ' 






\ 1 


4 












3 


{F'^BF'^B) 


2 


> ID 


1 


2 


> 


-8 










2 


-8 






{FB^FB"^) 


1 

X 

> 




1 


-1 

> 






x-' / 









' 












1 


-6 






{F^BF'^B) 


u 














{F^B'^) 


4 














{F^B^F^B^) 


2 




1 


2 




4 


{F^B^FB^) 


4 


► 35 








► -42 




(F^BF^B^) 


4 




2 


16 






{F^B^) 


2 




5 


-50 






(FB^FB^) 


2 




2 


-8 






(FB^FBF^B) 


4 

J 




1 


4 

> 





Table 2: Same as in the Table but using Eq. (|23|) 



The total number of nonvanishing terms on the RHS of Eq. (|2.16|) for the pure bosonic 
case, which we shall denote as is given by the following generating function 



n=0 



2n 



f Jl-4(m-l)c2-f + 1 



(2.39) 



1 — c^m^ 

This formula is derived in the Appendix A using non-commutative free random variables. 
For m = 2 Eq. (|]3|) yields 



-,2n 



1 



E 



(2n)! 



^2.40) 



These numbers describe the sum of the combinatorial numbers in the third column of the 
Tables | - 1. 



2.5 Relation to the Kazakov— Migdal model 

Equation ( ^.39 ) is known from the solution of the Kazakov-Migdal model with the Gaus- 
sian potential [|l^. There is the following reason for that. Suppose that the matrices Wa 
are unitary instead of the general complex ones. Then one has 

trC/„.£/i . . . £/„,._,C/t \ = j \ '-P^ . (2.41) 

\ iV ^" / Haar measure I for open loops 
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Structure 



{A^BA^B) 

(A^B^) 

{AB^AB^) 



{A^B^) 

(A^BA^B) 

{A^B^) 

{A^B^A^B^) 

{A^B^AB^} 

{A^BA^B^) 

{A^B"^) 

(AB^AB^) 

(AB^ABA^B) 



Combinatorics 



1 



1 
4 
2 
2 
1 



> 10 



6 
6 
4 
2 
4 
4 
2 
2 
4 



^35 



Value 



1 



2 
1 



5 
2 
1 
2 
1 



14 

5 

2 

4 

3 

2 

2 

5 

2 

1 



Contribution 



1 




25 

16 

2 

8 

1 



^52 



196 

150 

24 

64 

18 

16 

16 

50 

8 

4 



)■ 546 



Table 3: Same as in the Tables |i| and Q but using Eq. ( 2.16| ) at m 
of Wa are denoted as Wi = A,W2 = B. 



2. The components 



Here the loops represent the sequences of indices {oi, 02, ... , 02n-i5 a2n}. The nonvanish- 
ing result is only when the loop is closed and encloses a surface of the vanishing minimal 
area, i. e. each link of the loop is passed at least twice. This is analogous to the so-called 
local confinement in the Kazakov-Migdal model. 

The generating function (|2.39|) coincides with the following correlator in the Kazakov- 
Migdal model with the Gaussian potential on an infinite D-dimensional lattice: 



n=0 



2n 



(^tr0^(O) 



/ n. d<P{x) nfLi dU,{x) e -^[^'^1 i tr ^^(o) 



with the action 



D 



^2.42) 



(2.43) 



provided that 2D = m. The integration over the unitary matrices U^{x) goes over the 
Haar measure. 

The solution of the Kazakov-Migdal model with the Gaussian potential can be com- 
pletely reformulated as a combinatorial problem of summing over all closed loops of a 
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given length with all possible backtrackings (or foldings) included. Its solution is 
given by Eq. ( pISQl) . 

By virtue of the Eguchi-Kawai reduction []12|^; the correlator ( |2.42| ) in the Kazakov- 
Migdal model on the infinite lattice is equivalent to that in the reduced model given 
by 



n=0 



2n 



with m = 2D and the reduced action being 



S f/] = f tr 02 + I tr 0^ - ciV E tr (^if/l^sf/c 



a=l 



^2.44) 



(2.45) 



We have introduced here the index a running from 1 to m = 2D for the reduced model 
to distinguish from the index /x running from 1 to D for the Kazakov-Migdal model on 
the infinite lattice. The representation ( 2.44 ), (|2.45|) can be finally rewritten as 



n=l 



2n 



' / e Ea=i jj tr (jy,Ui(f>2U. 

/rf0irf(/)2e-^[<^'^] 



;2.46) 



Haar measure 



In order to prove the equivalence of Eqs. ( p.44| ) and ( p.46| ), we calculate the integral 
over dcpi and d(j)2 using Eq. ( [^.17] ) with Wa substituted by Ua- 



?2e 



det-i/2 



I®I 



c' E UaUl ® (u,ul) 



a, 6=1 



^2.47) 



The determinants in the numerator and the denominator on the RHS of Eq. (|2.46| ) ob- 
viously cancel before the averaging over UaS. An analogous cancellation happens in 
Eq. ( 2.44 ) as well in spite of the fact that each of them is averaged over its own UaS. The 
point is that the determinant ( |2.47|) under the sign of averaging over U behaves at large 
N [/-independent constant. 

The representation ( p.46|) looks very similar to the generating function ( p.l3| ) of the 
meander numbers. The difference is that the average is over the unitary matrices in 
Eq. ( p.46| ) and over the Gaussian complex matrices in Eq. ( p. 131 ). 

We can interpolate between the two cases by modifying the weight for averaging over 
W^s along the line of Ref. ^M. Let us introduce 



a=l 



w, 



(2.48) 



Then the averaging over the Gaussian complex matrices is reproduced as a ^ while the 
average over the unitary matrices is recovered as a ^ oo since the matrix Wa is forced 
to be unitary as a —>■ oo. 



'See Ref. |13|1 for a review. 
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We see, thus, that the words are the same both for the meander problem and for the 
Kazakov-Migdal model. The only difference resides in the meaning of nonvanishing words 
— it equals unity for the unitary matrices. 



3 Representation via non- commutative variables 

The set Ua, u\ of non-commutative variables obey the Cuntz algebra 

Uaul = Sab ■ (3.1) 

It is convenient to consider them, respectively, as annihilation and creation operators in 
a Hilbert space with the vacuum \Q) which satisfies 

= 1. (3.2) 

The completeness condition says that 

m 

Y^uiua = i-\n){n\ . (3.3) 

a=l 

There are no more relations between the non-commutative variables. 
3.1 Bosonic case 

Let us construct the generating function for words via the non-commutative sources Ua 
as 

GM = 1 4tr ^ 



°° 1 /I \ 

= E Ua,Ua,---Ua,(—tlWa,Wa2---Wa,,J (3.4) 

where the average over a = 1, . . . , m matrices Wa is with the Gaussian weight as before. 
Then the generating function (|2.13| ) for the meander numbers is given by 



9.\Gx{u)Gx{u^)\^) = c + cmM{c) (3.5) 



with c = 1/A^. The contraction of indices on the RHS of Eq. (|2.3|) is obviously reproduced 
using Eqs. (|l|) and (|3). 

The generating function Gx{u) obeys the Schwinger-Dyson equation 

m 

XGx{u) - 1 = E Gx{u)uaGx{u)ua , (3.6) 



a=l 
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which can be derived in a usual way, by shifting W. The cychc symmetry of the trace 
imphes 



XGxiu) -l = J2ub (XGxiu) -l)ul (3.7) 

b=l 

Inserting here Eq. (|3.6|), we get 

m m 

J2 Gx{u)UaGx{u)Ua = UaG x{u)UaG x{u) . (3.8) 



a=l a=l 



An alternative combinatorial derivation of Eqs. ( p.5| ) and ( |3.6| ) is presented in the Ap- 
pendix B. 

We shall use for Gx{u) and Gx{u^) the short-hand notation 

Gx = Gxiu) , G{ = Gxiu^) , (3.9) 
so that Eq. ( p. 61) can be rewritten as 

XGx - 1 = GxUaGxUa , 

XGi-1 = ulGWaGi, (3.10) 

where the summation over repeated indices is implied here and below except when it is 
specially indicated. Using Eq. (p.l|), one alternatively rewrites Eq. ( |3.10D as 

(XGx -l)ul = GxUaGx , 

uJXGi~l) = GWaGl (3.11) 



If m's were ordinary commutative variables, the solution to the quadratic equation 
would be simple 



2 u 

For 171 = 1 this is nothing but Wigner's semicircle law and one reproduces Eq. ( ^.18] ) by 
expanding in 1/A. 

A formal solution for the non-commutative variables can be obtained representing 



Eq. ( |3.6|) as 

Gx = T (3.13) 



Iterations of this equation, as was found by Cvitanovic ||15[, lead in the continued fraction 



Gx{u) = ^ . (3.14) 

A - Ua^ Y 



X '^02 ]^ 
A ~ 
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Expanding the RHS in 1/A^, one gets 

GaN = ^ + ^ + 



+ + 0[-^) (3.15) 

where stands for UaUa- All possible planar combinations of n's appear to next orders 
(with unit coefficients) while the total number of terms to the order A~^"~^ equals the 
Catalan number Cn- Technically, it is more simple to derive the expansion of G\ (u) in 
1/A by direct iterations of Eq. ( |3.6|) . The substitution into Eq. ([Ol ) then recovers the 
lower meander numbers. 



Equation (p.6| ) is in fact well-known for the Gaussian models. It is a consequence 
of the relation between the generating functionals for all planar graphs Gx {u) and for 
connected planar graphs 



oo m , 

E 3aJa,---3a,Aj;^^^Wa.Wa,---Wa,„) (3.16) 

n=0 ai,a2,---,a2n=l 



which says [|l^, [T^ 

a (u) = W (XuG (u)) (3.17) 

while the cyclic symmetry gives 

W {XuG (u)) = W {XG (u) u) . (3.18) 
There is only one connected graph in the Gaussian case so that W (j) is quadratic 

W(i) = ^ + ^. (3.19) 

Equations ( p.l7|) and (|3.18|) now recover Eqs. (|3.6| ) and (|3.8|) , respectively. 

It would be interesting to apply the theory of non-commutative free random vari- 
ables [0, whose application to matrix models has been discussed recently in Refs. |1^, |T^, 
for this problem. 



3.2 General case 

For a general complex matrix model when some components of Wa are bosonic and some 
are fermionic, we define the generating functions 

Gx(u) = /iti- ^ 



N X^-EZb=l^aUtWaW,/^^^^^ 

oc ^ m , ^ _ _ . 

„=0 A'"+\i,a2,-,a2n=i /Gauss 
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and 

CO m I ^ _ _ > 

E E ^«i^«2 • ■ ■ 1*a2„_iW,,„ ( — tr ■ ■ ■ l^a2„_iW^a2„ )^ (3.21) 

n=0 ^ ai,a2,-,a2„=l ^ ^ GauSS 

where the non-commutative sources are the same as before. 

The generating functions G\{u) and G\{u) are not independent due to Eq. ( p.32|) . 
One gets 

GaK) =GA(ixJ (3.22) 

where 



Ua = ^o{a)Ua (3.23) 

and a(a) is defined by Eq. ( p.30 ). In other words the components of u^, which are associ- 
ated with fermionic components of Wa should be multiphed by i on the RHS of Eq. (|3.22|) 
while the ones for the bosonic components remain unchanged. 

The Schwinger-Dyson equations for G\{u) and G\{u) can be derived in a usual way 
from the recurrence relations 



jV / Gauss 

= E'^'^2„a2...(4t^W^«i^«.---W^«2.)^ (h^'^a,,,,■■■W^,^^ (3.24) 



fc=0 

and 



-iV / Gauss \ -iV / Gauss 



1 



ESa.a,Al^^'Wa.---Wa,.-y ( W^,.,. ' ' ' Wa,.) ^, ■ (3-25) 



k=l 



N "^'=-VGauss\iV ™ "^V Gauss 



They can be obtained in a usual way by shifting W2n and Wi, respectively. 

Equations (|3.24|) and (|3.25|) and the cyclic symmetry of the trace result in the equations 



XGx{u) - 1 = E UaGx{u)uaGx{u) = J2 Gxiu)uaGxiu)ua , (3.26) 

a=l a=l 

and 

m m 

XGxiu) - 1 = E (^ia)Gxiu)uaGxiu)ua = J2 (y{.ci)uaGx{u)uaGx{u) . (3.27) 

a=l a=l 

These two equations are not independent due to the relation (p.22|) . One can be obtained 
from another by the substitution 

Ua^Ua, (3.28) 
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where Ua is given by Eq. (|3.23| ). Their combinatorial interpretation is presented in the 
Appendix B. 

Introducing in addition to Eq. ( |3.9| ) the short-hand notations 

= Gxiu) , G{ = G\iu^) , (3.29) 
one derives from Eqs. ( p.26[ ), ( |3.27| ) the analog of Eq. ( |3.11| ) 



as 



{\Gx-l)ul = G\UaGx. 



(3.30) 



and 



\Gx - l) 1*1 = a{a)GxUaGx , 
= a{a)G\uiG{. 



a \^G\ — 1 



(3.31) 



The solution to Eqs. ( p.26| ) and ( p.27| ) can be easily obtained for a pure fermionic 
model and commutative sources. Summing up Eqs. (|3.26|) and (|3.27]), we get 



Gx — ^ — Gx 



\Gx-l= u^G\ - ^Gx 



2 ^ 

a' 



(3.32) 
(3.33) 
(3.34) 

2 ■ .,.2V'> . ■ (3.35) 

Equations ( ^.331 ) and ( |3.34| ) for m = 1 recover the ones of Ref. while the expansion of 
Eq. i ^I^ in 1/A gives Eq. ( ^lOl) . 

In order to find a formal solution to Eqs. ( |3.26|) and (|3.27|) in the form of a continued 
fraction, let us rewrite them as 

1 



and 

with the solution 
and 



Gx = \ + :^--^^X' + Au' 



and 



G), 



X - UaGxUa 

1 



(3.36) 



(3.37) 



A - UaGxUa 

where Ua is defined by Eq. (|3.23| ). Iterations of this equations lead to the following analog 

of ( CT) 

Gx [u) = ^ (3.38) 



A - Wai ■ 



\ — u 



1 



a2 



0.2 



X - u„ 



1 "^(13 
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and ^ 

Gxiu) = ^ . (3.39) 



A - Wai ■ 



A - - 



A '^03 "^(13 



Here u and w interchange in the consequent hnes of the continued fractions. 
The quantities on the LHS's of Eqs. (|3.3CI|) and (|3.31|) read exphcitly 



00 -| m /I 

(AGl - 1) = E E <■■■ ( IT^tr W^W., . . . W^,^_,W^,, 

n=l a2, -,a2„=l ^ ^ GauSS 

(3.40) 

and 

(AGa - 1) = E E ^*ailia2---Wa2„_i ( ^triy,,iy„, 

^ ^ n=lA^ ai,a2,.,a2„_i=l / GauSS 

(3.41) 

The generating function (|2.2tj|) is determined by 



n\Gx{u)Gx{u^)\n) = {n\Gx{u)Gx{u^)\n) = c + c{mb - mf)M{c) (3.42) 



where mh and rrif are the numbers of bosonic and fermionic components of Wa, respec- 
tively, and A^ = 1/c. This formula follows from the fact that the LHS reproduces the 
contraction of indices as in Eq. (|2.25|) using Eqs. ( p.lj) and ( |3.2|) . 



3.3 Supersymmetric case 

For the supersymmetric case m?, = ruf = 1 and Eq. ( p. 421) does not determine the meander 
numbers. One should use, instead, Eq. (|2.37|) or Eq. ( p.38|) , where there is no summation 



over one of the indices, to get the principle meander numbers. 
Let us denote the components of Ua as 

Ua = {u,v) . (3.43) 

Equations ( |3.40| ) and ( |3.41| ) then result in 



^A 

or alternatively 



M(^) = (Q\Gxu^uG{\Q) = -A^ (n\Gxv^vG{\n) (3.44) 



M(^) = A^ (n\Gxu^uG{\n) = -A^ (n\Gxv^vGi\n) . (3.45) 
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The equality sign between the two expressions on the RHS of Eq. (|3.44| ) or Eq. ( |3.45|) 
is due to the supersymmetry. This can be shown using the completeness condition ( |3.3| ) 
which takes in the supersymmetric case the form 

u^u + v^v = 1-\Q) {Q\ . (3.46) 

Inserting this, say, between G\ and G\, we get 

X^Gxu^uG{ = -X^Gxv^vG{ + X^GxG{ - \Q) {Q\ . (3.47) 

Taking the vacuum expectation value of this formula and remembering that 

{Q\ GxG{ \Q) = 1 (3.48) 

due to the supersymmetry, we prove the statement. 

Using Eqs. (p:30|) and (g), Eqs. j ^J^ ) and (|^) can be rewritten as 

M(^) = (n\GxuGxGWGi\n) = (n\GxvGxGWG{\n) (3.49) 

and 

M(— ) = (n\GxuGxG{u^G{\n) = (n\GxvGxG{v^G{\n) . (3.50) 



•A2- 

This representation of the generating function for the principle meander numbers is con- 
venient for an iterative procedure. 

To perform calculations order by order in c = it looks reasonable to try to use 

the fact that there are only two non-commutative variables and to expand in v. This is 
a standard trick of dealing with two non-commutative variables in the Hermitean two- 



matrix model 20, 21 



We introduce, therefore, the expansions 

oo 

Gx = EG*^"^ 



A 5 

n=0 

oo 

in) 



Gx = E(-irGA , (3.51) 



n=0 



where (—1)" in the expansion of Gx is due to Eqs. ( |3.22| ) ( |3.23| ). We have 

Gl"' - '-^Zr^ (3.52) 



while G^x^ involves v exactly 2n times. As we shall see below, 

gI"^ ~ TT-TT for ^ odd , 

^ \2n+l ' 



gI"^ ~ forn>2even. (3.53) 
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Therefore, a contribution to M„ arises in Eq. (|3.45|) at most from G^"^ for odd n and from 
G^x even n. 

The functions G^^^ can be found recursively from Eq. (|3.36| ) which takes the form 

(3.54) 



A ^ 



1 - E:f=i(-l)"«Gt\Gf - E^=oi-^rvGi'^vGf^ 
At each step one has to solve the equation 

^(n) 



(3.55) 



with some whose solution is given by 



oo , 



(0) 



u 



1=0 

oo 



Gf> = i:(Gl%)'^,,(Gf.)', 

1=0 



(3.56) 



where p > 1. Both formulas can be rewritten in a unique way via the contour integral 



G 



(n) 



duj 



1 



2m 



CJ- (-l)"Gf M 



1 



1 - uGf^u 



(3.57) 



We shall also introduce the obvious short-hand notation for RHS's of Eqs. ( |3.56| ) or ( |3.57|) 
as the brackets 

4"^ = {An} . (3.58) 



The functions An can be found recursively from Eq. ( [3.541 ). Few lower ones read 
explicitly 



Ai 

A2 

A. 



GfKGf^Gf^ , 



1 



G^. 



(1) 

A ? 



G^'^G^'^vG^'^ + G 



(2) 



:G 



(1) 



G 



(1) 



-G 



(2) 



'-'A ^(0)"-^A (0)' 



G 



G« . (3.59) 



They obviously satisfy the equation 



n— 1 



XAnV^ = J2i~^)'G 

k=0 



fc^(n-l-fc)^^{fc) 
A A 



and analogously 



n-l 



^vAi = Y.{-lfG\^'^v^G\^-- 



l-k) 



(3.60) 



(3.61) 



fe=0 



which can be obtained substituting the expansion (|3.51| ) into Eq. ( |3.30|) (or Eq. ( p.31| )) 
and using the fact that 

(3.62) 



G^^V^ = AnV^ 
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due to Eq. (|336|) . 

Using Eq. (|338| ), we get from Eq. (|339|) for gJ^'^'s 



= -Gf {t^Gf ^Gf + Gf^{{vG^;:^Gf\vGf\:}}Gf^ , 

Gf = -Gf {t;Gf {^Gf {^Gf4^A°^^}G'f^^}Gf 
+Gf{.Gr{.Gi%}Gr{.Gl%}Gl%}Gf 
-Gr{{.Gi%}Gf{.Gf{.Gi%}Gi%}}Gr 
+Gf{{.Gi%}Gf{.Gi%}Gr{.Gi%}}Gf 

-Gf {{.Gf {.Gi%}Gi%}Gr{.Gi%}}Gr , (3.63) 

where we have used the fact that G^^ commutes with a bracket but not with v. There 
are no problems here with the order of brackets which act hke in TeX: they immediately 
contract each other when it is possible. This contraction is associated with the summa- 
tion (|3.56|) (or the integration (|3.57|) ). The general G^"'' contains G„ (the Catalan number 
given by Eq. (|2.5| )) terms of this kind with alternating signs. Some rules for representing a 
general term can be formulated which resemble the Wick pairing of bilinear combinations 
of f 's. 

It is worth noting that most of the above formulas would look quite similar for the case 
of just 2 bosonic fields, i.e. of m = 2, while there will be no alternative signs for 2 bosons. 
Due to these minus signs, the supersymmetric case is somewhat simpler than the m = 2 
one. Say, the leading-order terms are cancelled in Eq. ( |3.59| ) for A2 and in Eq. ( 3.63 ) 



(2) 

for G'^'. This is a reflection of the general property ( p^.l9D of Gaussian averages of the 



fermionic matrices. It is the reason why the leading order term vanishes in Eq. ( p.53| ) for 
even n. In addition some relations are imposed by the supersymmetry. The simplest one 
follows from the substitution of the expansion ( p.51|) into Eq. ( p.48|) which yields 



00 

A' E(-l)" (^1 ^1"^^!^"^ \n) = l. (3.64) 

n=0 

Only the diagonal terms survive here due to the definitions (|3.1[), ( |3.2| ). 

In order to calculate M„, we substitute the expansion (|3.51| ) into Eq. ( p.45| ) which 
gives 

-1 00 

^(T^) = A'E(-l)^""'M^l4"Vt;Gi(")|n) . (3.65) 



^ n=l 



We can also rewrite the RHS as 



n=l 

due to Eq. (|33^) . 



-| 00 

Miy,) = E(-l)^""'^ (^1 ^nv^^^l \n) (3.66) 
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It is easy to calculate the term with n = 1 in Eq. (|3.65| ) (or Eq. ( |3.66|) ). Using 
Eq. ( p.63| ), we get 

G'i^v^ = ^Gf , vGl^''> = jG{^'^ gI'^'^ , (3.67) 
so that Eqs. (^), (U) and result in 

M{y2) = ^HT.c''CI] +A^5:(-l)("-i)(fi|Gl"V.Gi(")|fi) . (3.68) 

^ \k=0 / n=2 

The first term on the RHS recovers Mi = 1 and M2 = 2 while the contribution of the 
next terms are controlled by Eq. ( ^.53 ). 

Let us demonstrate of how this iterative procedure works by the explicit calculation 
up to = 3 when at most G^p is essential. We get from Eq. ( p.63|) 



(2) -j- _ —vuv'^u + uv'^uv f ^ \ ^t(2) _ v^u'^v'^'^u'^ — u^v'^'^u^v'^ n ( ^ 

GlV^4 + 0(i). „Gr.-!^ + o(l). (3.69) 

The substitution into Eq. ( |3.68[ ) recovers M3 = 8. 

A most difficult part of the described general iterative procedure is to calculate 
{Vt \ G^'^^'u'^f G^^''"^ \VL) which involves terms for n > 1. There is no problem to calculate 
a contribution from each individual term by using the formula 

{Vt\ fQ{u)vfi{u)v ■ ■ ■vfn{u)gn{u^)v^ ■ ■ ■ v'^ gi{v))v'^ gQ{u^) \Vt) 

= Uu)g^{v)) |^]) h{u)gr{v)) |^]) • ■ ■ Uu)gn{u^) |^]) , (3.70) 



which is a consequence of Eqs. ( P . 1| ) and ( |3.2| ). When calculating order by order in 
c = it is more convenient first to collect similar terms in G^^^f^ and fG^*'"'' to 

the given order c". After this each term in G^"^f^ is orthogonal to all terms in vG^^"'' 
except for one. Whether or not this kind of an iterative procedure could be useful for 
deriving recurrence relation between meander numbers would depend on a possibility to 
reexpand (|3.59|) in such orthogonal terms in the general case. 

We see now the difference between the expansion in v for the meander problem, de- 
scribed in this Subsection, and that for the two-matrix model with a polynomial poten- 



tial pO[. The latter reduces to an algebraic equation while the former does not. This 



differs the matrix models describing the meander numbers from the Kazakov-Migdal 



model whose solution can be described in the language of the two-matrix model |22 



4 Discussion 

We have considered in this paper the meander problem which results in more complicated 
matrix models than those solved before. It belongs to the same generic class of problems 
of words as, say, the large- QCD in = 4 but is presumably simpler. 
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We have introduced a new representation of the meander numbers via the general 
complex matrix model. This allowed us to consider the supersymmetric matrix model 
which includes both bosonic and fcrmionic matrices and describes the principal meanders. 
This model looks simpler than a pure bosonic model and one can try to solve it by using 

supcrsymmctry Ward identities. 

Using non-commutative sources, we have reformulated the meander problem as that of 
averaging in a Boltzmannian Fock space whose annihilation and creation operators obey 
the Cuntz algebra and have shown the equivalence with the combinatorial approach based 
on the arch statistics. The averaging expression is represented in the form of a product 
of two continued fractions but further progress requires their better understanding as 
functions of the non-commutative variables. One of the possible approaches can be based 
on the standard procedure of disentangling via the path integral. 

We have discussed also the relation between the matrix models describing the meander 
problem and the Kazakov-Migdal model on a D-dimcnsional lattice. The words are the 
same both for the meander problem and for the Kazakov-Migdal model while the only 
difference resides in the meaning of nonvanishing words. This relation could give a hint 
on how to solve the meander problem. 

We have demonstrated how the solution of the Kazakov-Migdal model with the 
quadratic potential can be obtained using the theorem of addition of free random vari- 
ables. We have shown that this approach does not work for the meander problem even 
within the interaction representation since the variables are not free for this case so that 
the theorem of addition is not apphcable. This observation casts doubts that a solution 
of large- A?" QCD can be obtained in the language of a masterfield given by free random 
variables. 

The supersymmetric matrix models of the type discussed in this paper for the meander 
problem are novel ones. It is worth studying them in connection with other physical ap- 
plications, in particular, with a discretization of super-Riemann surfaces and superstrings. 
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Appendix A Derivation of Eq. ( |2.39| ) via free ran- 
dom variables 



Equation ( |2.39|) can be derived using free random variables as follows. First of all let us 
note that 



#n = # of nonvanishing terms in ( ^ tr I ^ W^W^J 

^ \a,b=l 



Gauss 



2n 



# of nonvanishing terms in ( ^ ^r | Aj 



\i=l 



(A.l) 



where the matrices Ai are Hermitean. The simplest nontrivial example is that of m = 
2, when the LHS has for n = 2 six nonvanishing terms of the type of WiWlWiWi , 
WiW^W2Wl W1WIW2WI W^WlW^Wl, W2WIW1WI W2W^WiWl The RHS has six 
nonvanishing terms of the type of A"^, A^B^, AB^A, B^, B^A^, BA^B, where we have 
used the notations Ai = A, A2 = B — the same as in the Table 3. 

If we now introduce free random variables A,-, such that 



n 



A 



2n 



A 



2n+l 



1, 

0, 



then 



2n 



\i=i 



n 



(A.2) 



(A.3) 



since the vacuum average of each nonvanishing term equals one. The resolvent for each 
A,- reads 



n 



A^ 



n 



n=0 



yU+l 



E 

n=0 



1 



2n+l 



(A.4) 



The inverse function to (|A.4|) is given by 

A{z)^Rr\z) = 



VI + 4^2 + 1 
2~z 



(A.5) 



Then the master-field operator can explicitly be constructed as 0, [19 

./l + 4(at)2 - 1 



2at 



(A.6) 



where the RHS is understood as the expansion around zero. 

The theorem of addition of free random variables states that if Aj's are free random 
variables and so they are representable as 

a, + Mal), (A.7) 



A, 
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where /j's are functions analytic around zero, then 



(A.8) 



Using this theorem, we get for the m matrices 



1 



A{z) = - + m[Ai{z) - - 



m 



VI +4^2 + 2 



m 



2z 



(A.9) 



The resolvent is given by the inverse function 

(2 — m)z + z^ — 4(m — 1) 



R{z) = A-\z) 



2 (z^ — rn? 



(A.IO) 



Equation (|2.39|) can be reproduced substituting z = l/c and dividing by c. 

As far as the Kazakov-Migdal model is concerned, Eq. (|2.46|) results in the following 
representation of the correlator (|2.42|) via the words built up of the m = 2D unitary 
matrices: 



-tr02(O)) = E C 



2n 



2D 

E 

n=0 a,i,a2,---,a2n- 



I 1 



Haar measure 



.(A.ll) 



Since the meaning of each word is either one or zero, the square on the RHS is not 
essential. Equation ( |2.39| ) with m = 2D then yields 



/ ^tr02(o) 



D^1-A{2D- I)c2-D + 1 
1 - 4L)2c2 



(A.12) 



which reproduces the solution of Ref. 10 



We discuss in the Appendix C that the approach, which is described here for the 
Kazakov-Migdal model, does not work for the meander problem since the random vari- 
ables are not free in the latter case. 



Appendix B Combinatorial interpretation of the re- 
sults of Sect. 3 

The representation of Sect. 3 of the meanders via non-commutative random variables 
can be alternatively derived pure combinatorially which clarifies the relation between our 
approach and that of Ref. . 

Equation ( p.6|) can be derived combinatorially as follows. Let us denote by Qn the 
sum of all possible arch configurations of order n with m colorings. For example, Qi for 
m = 2 is given by 
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Qi = ( \ + f \ , (B.l) 

where the two terms on the RHS differ by the coloring. Picking up the leftmost arch, we 
have the following recurrence relation for QnS: 

^" = EE ^n-k-l, Qo = ^- (B.2) 

fi^=Q colors 

of n 

Let us analogously denote by the sum of all possible upside-down arch configura- 
tions of order n with m colorings. For example, QI for m = 2 is given by 

Ql obeys the recurrence relation which is similar to Eq. (B]2). 



Let us now introduce the "multiplication" o of arch configurations. Its meaning is 
evident from the following examples: 



We have then the relation 



J the sum of all possible 1 
I multi-component meanders J 



Sn°Gl= < , } , (B.4) 



and 

n 

# of terms in ^„ o = ^ M^^^m'' , (B.5) 

fc=i 

since each loop in the meander can be colored with m colors. 

The "multiplication" o can naturally be represented in the operator language as fol- 
lows. Let us introduce m non-commuting variables Ua (a = 1, . . . , m), one for each color 
in the arch configurations. Let us associate with each arch configuration a word made up 
of in a way, which is evident from the following example: 

< — > U1U1U2U2, if ui ^ — and U2 — . 

Analogously, we associate words, made up of xt^'s, with the upside-down arch configura- 
tions as in the example 
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(niitiW2it2)^ = U2U2u\u\, if Ui ^ — and U2 ^ — 



We have from Eq. (p.2|) the following equation in terms of UaS: 

n—l m 

Gn{u) = UaGk{u)UaGn-k-l{u) , 

k=0 a=l 



(B.6) 



where G„('u) stands for the sum of the words of order n. Introducing the generating 
function (cf. Q) 

00 

Gxiu) = J2Gniu)X~'--\ (B.7) 



n=0 



Eq. ( p.6| ) can be rewritten in the form ( p.6|) . 
There is the correspondence 



AoB 



(B.8) 



Using this correspondence, Eqs. ( |B.5| ) and ( |B.7D , we arrive at Eq. 

Equation ( |3.26| ) can be interpreted in the same way. As an example, let us consider 
the supersymmetric case of m = 2, when Ua = {u,v). The analog of Eq. ( [B.2| ) for this 
case reads 



Gn 



n—l 



[ ( Gk \ Gn-k-l + ( Qk \ Gn-k- 



k=0 



(B.9) 



where 



Gn = Gn{n,n), 
Gn = Gnin,-r\) 



[B.IO) 



Let us introduce the projector V, such that it picks up the arch configurations whose 
rightmost arch is of light color. For example. 



Then, we have 



+ I r^ 



the sum of all principle 
meanders of order n 



V[Gn]oV Gi 



(B.ll) 



The projector V corresponds in the operator language to the insertion of m^u in the 
Eq. (1^. 
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Appendix C Comment on freeness of random vari- 
ables for the meander problem 

We demonstrated in the Appendix A how the Kazakov-Migdal model can be solved via 
free random variables. We comment here that this approach does not work for the meander 
problem since the random variables are not free for this case. 

The meander numbers can be represented by the vacuum average of certain non- 
commuting variables as 



Gauss 



\i=l 



(C.l) 



where the operators Aj's and Aj's under the vacuum averaging are 



Ai — Cli ~\~ Oij 



A. 



ai + a] . 



Here the two sets of operators a^'s and Oj's obey the Cuntz's algebra 



t _ r ~ ~t _ 



(C.2) 



(C.3) 



independently of each other while all the operators with tildes commute with all the 
operators without tildes. The equivalence of the two representations in Eq. ( |C.1| ) is 
evident from the factorization 



(Tr(A® B)) = (tr AtrS) = {tr A) {ti B) 



(C.4) 



at large A^. 

Equation (|C.1| ) involves the direct product of matrices which is represented as the 
product of two independent operators Ai and A^. The non-commuting variables of the 
type AiAi are not free random variables of Ref. [1^ since they do not satisfy the defining 
axiom of free random variables. This can be seen by considering the average 



{AAf - {^{AAy 



- ({B^y 



{BBf-{{BB)' 







(C.5) 



which would vanish if they were free random variables. 

It is easy to see by direct calculation that the expression ( |C.5| ) does not vanish. We 



rewrite it as 



n 



{AAy - 1 



:bb) 



(AAf 



(BBf 



n 



A' 



n 



n 



B' 



n 



n) +1 



3^-2^-2^ + 1 



(C.6) 
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since 



n 



(AAf 



n 



n 



{BBf 



n 



1. 



(C.7) 



The squares on the RHS of Eq. (|C.6|) are due to analogous contributions from the operators 
with tildes. 



The nonvanishing RHS of Eq. ( |C.6| ) means that the variables AiAi are not free. This 
is, as is already mentioned, a consequence of the direct product of matrices which results 
in the product of Aj and Ai and leads, in turn, to the squares in Eq. ( C.6 ). If there 
are no squares, the axiom of freeness is obviously satisfied. The same is true for the 
supersymmetric well. 

The example of this Appendix shows that random variables are not necessarily free 
in (multi-) matrix models with complicated interaction even within the framework of the 
interaction representation. The simplest average for the meander problem which violates 
the axiom of freeness is that ( |C.5| ). For this reason, the theorem of addition with 



n=0 



{AiAi 



E 

n=0 



2n+l 



(C. 



where C„ are given by Eq. (^.5|), can be used, quite similarly to the Appendix A, to 
describe the meander numbers only up to n = 3 but fails to reproduce M„ for n > 4. 



30 



References 



[1] J. Ambj0m, B. Durhuus and J. Fr51ich, Nucl. Phys. B257[FS14] (1985) 433; 

F. David, Nucl. Phys. B257[FS14] (1985) 45; 
V.A. Kazakov, Phys. Lett. 150B (1985) 282; 

V.A. Kazakov, I.K. Rostov and A.A. Migdal, Phys. Lett. 157B (1985) 295. 
[2] E. Marinari and G. Parisi, Phys. Lett. B240 (1990) 375. 

[3] L. Alvarez-Gaume and J.L. Manes, Mod. Phys. Lett. A5 (1991) 2039; 

G. Gilbert and M. Perry, Nucl. Phys. B364 (1991) 734. 

[4] V.A. Kazakov and A.A. Migdal, Nucl. Phys. B397 (1993) 214. 
[5] V.A. Kazakov and I.K. Kostov, unpublished. 

[6] P. Di Francesco, O. Golinelli and E. Guitter, Meander, folding and arch statistics, Saclay 
preprint SPhT/95-059, |hep-th/9506030 . 

[7] S. Lando and A. Zvonkin, Theor. Comp. Sci. 117 (1993) 227. 

[8] Yu. Makeenko and K. Zarembo, Nucl. Phys. B422 (1994) 237. 

[9] G.W. Semenoff and R.J. Szabo, Int. J. Mod. Phys. A12 (1997) 2135. 

[10] D. Gross, Phys. Lett. 293B (1992) 181. 

[11] Yu. Makeenko, Mod. Phys. Lett. A8 (1993) 209. 

[12] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 1063. 

[13] S.R. Das, Rev. Mod. Phys. 59 (1987) 235. 

[14] E.-M. Ilgenfritz, Yu.M. Makeenko and T.V. Shahbazyan, Phys. Lett. 172B (1986) 81. 

[15] P. Cvitanovic, Phys. Lett. 99B (1981) 49. 

[16] E. Brezin, C. Itzykson, G. Parisi and J.B. Zuber, Commun. Math. Phys. 59 (1978) 35. 

[17] D.V. Voiculescu, K.J. Dykema and A. Nica, Free Random Variables, AMS, Providence 
1992. 

[18] M.R. Douglas, Phys. Lett. B344 (1995) 117. 

[19] R. Gopakumar and D.J. Gross, Nucl. Phys. B451 (1995) 379. 

[20] M.R. Douglas and M. Li, Phys. Lett. B348 (1995) 360. 

[21] S.M. Carroll, M.E. Ortiz and W. Taylor IV, Nucl. Phys. B468 (1996) 383. 

[22] M.I. DobroUubov, Yu. Makeenko and G.W. Semenoff, Mod. Phys. Lett. A8 (1993) 2387. 



31 



